Increasing d-wave superconductivity by on site repulsion 
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We study by Variational Monte Carlo an extended Hubbard model away from half filled band 
density which contains two competing nearest-neighbor interactions: a superexchange J favoring 
d-wave superconductivity and a repulsion V opposing against it. We find that the on-site repulsion 
U effectively enhances the strength of J meanwhile suppressing that of V , thus favoring supercon- 
ductivity. This result shows that attractions which do not involve charge fluctuations are very well 
equipped against strong electron-electron repulsion so much to get advantage from it. 
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The interplay between strong correlation and super- 
conductivity is one of the major problems raised by the 
discovery of cuprate High-Tc materials. 

Indeed, within the conventional BCS theory for 
phonon-mediated weak coupling superconductivity, 
a strong electron-electron short range repulsion, 
parametrized by an Hubbard [/, can only depress Tc. 
Landau Fermi liquid theory identifies two main sources 
for this reduction; namely, as U increases, the quasipar- 
ticle wave function renormalization Z < 1 diminishes, 
meanwhile the effective mass increases. Thus each 
interaction amplitude, including the phonon-mediated 
attraction g, acquires a renormalization factor times 
a vertex renormalization. If the latter is negligible, the 
bare amplitude is reduced to f/, = Z^g. Therefore, as U 
increases, the dimensionless coupling which parametrizes 
the bare attraction, Ao — po\g\, being po the uncorrelated 
density of states at the Fermi level, is renormalized into 



A, — I (7* 



Z^ — Aq < Ao, 
m 



being usually Z < m/m^. On the other hand the 
Coulomb pseudopotential fi^ increases, so that A, — /i, 
diminishes even more pushing Tc down. 

By solving a model for alkali doped fuUerenes within 
dynamical mean field theory, it has been recently 
argued Ij that there exist attractive channels for which 
vertex corrections may compensate the wave-function 
renormalization factor leading to 

Ao Ao > Ao, 
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which may indeed lead to an enhancement of T^ by in- 
creasing U. The main ingredient was identified into a 
pairing mechanism not involving the charge density op- 
erator, which is mostly subject to the renormalization 
induced by U, but other internal degrees of freedom, like 
the spin (or the orbital index, if orbital degeneracy is 



present) , unveiling a kind of spin-charge separation even 
within Landau Fermi liquid theory. 

This proposal is not far in spirit from the original Res- 
onating Valence Bond (RVB) scenario for High-Tc su- 
perconductivity in the t-J model There supercon- 
ductivity occurs naturally upon doping since the par- 
ent insulating state is well described by an RVB state: 
the spin-singiet valence-bond pairs naturally evolve into 
Cooper pairs. They can propagate around the lattice 
only through the empty sites left behind by the holes. 
This constraint easily explains a superfluid density pro- 
portional to the hole doping. Moreover, although at 
small doping superconductivity is suppressed, the energy 
scale associated to the binding energy of the valence- 
bond pairs remains finite, which is advocated to explain 
the experimentally observed pseudogap phase of High-Tc 
materials(3, ^^l. 

Within the Fermi liquid framework, the constraint of 
no double occupancy appears to renormalize the quasi- 
particle hopping to a value Zt ~ 2St, S being the doping, 
while leaving untouched the quasiparticle attraction, here 
provided by the superexchange J. The superconducting 
phase of the t-J model can be approached either from the 
half-filled antiferromagnetic Mott insulator upon increas- 
ing doping or at finite doping by decreasing temperature. 
In both cases, even though the T = superconductor 
might still be described in terms of Landau-Bogoliubov 
quasiparticles, in the RVB language spinon-holon con- 
fined objects, the relevant excitations above Tc or in the 
close vicinity to the half-filled antiferromagnet do not 
necessarily look as conventional quasiparticles. 

For this reason, in this work we shall try to under- 
stand whether this strongly correlated d-wave supercon- 
ductor can be approached at zero temperature starting 
from a weakly correlated regime, where standard many- 
body techniques and the well established Landau Fermi 
liquid theory should apply. 

We consider an extended Hubbard model in two- 
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dimensions for the average number of electrons per site 
1 — (5 < 1 , namely away from half-filling, 

H = -t E E ('^Is- + ^•^•) +t^E"*T^a 

<ij> a i 

+j E ■ ^ + ^ E ^'^3' 

<ij> ^ ^ <U> 

where, in addition to the on-site repulsion, we add a 
nearest neighbor spin-exchange and a charge-repulsion 
term, with strengths J and V, respectively. These 
nearest-neighbor interactions compete, J favoring a d- 
wave singlet pairing away from half-filling while V op- 
posing against it. For V = Q and U strictly equal to 
oo, 1^ reduces to the standard t-J model, which also 
corresponds to the large U limit of the pure Hubbard 
model, in which case J — > At^ /U. However, contrary 
to the latter, model for J > y is undoubtedly a d- 
wave superconductor at weak coupling (U , V and J all 
much smaller than t) also within the Hartree-Fock ap- 
proximation. For this reason model ^ is more suitable 
to address the issue of electron-electron correlation effects 
on d-wave superconductivity. Moreover, since V involves 
charge-density while J spin-density operators, the pres- 
ence of both gives the opportunity to test if U indeed 
induces different renormalization factors on charge with 
respect to spin vertices. 

A variational approach which was shown to correctly 
reproduce both the weak |^ and the strong cou- 
pling limits of the Hubbard model appears well suited 
for model too. It consists in searching by a varia- 
tional Monte Carlo (VMC) technique for the best wave 
function of the form 

I*) = APN^^Jastrow^^Gl^SCs) (2) 

where A is a normalization factor, \^ bcs) ^ BCS wave- 
function projected by Pn onto a fixed number of par- 
ticles, with a d-wave gap- function = Ayaricosk^: — 
cosfcy), being A^^^ a variational parameter. Pq is a 
Gutzwiller projector: 

-Pg = ]^(1 -ao"'n,T"n,i)' (3) 

n 

whereas Pjastrow a long-range Jastrow factor which en- 
forces the correct long- wavelength behavior of the density 
structure factor: 

^Jastrow = e ^<'^> ^«^J» , (4) 

where ". . ." stands for the summation over next, next- 
next, etc. nearest neighbor sites (i.e. all those possible 
on a finite size sample). 

The method is based on the Stochastic Reconfiguration 
(SR) technique 6], which allows to minimize the energy 



of a variational wave function containing even a large 
number of parameters. 

To get further insight from the numerics, we compare 
the results with those obtained by the Gutzwiller Ap- 
proximation (GA) for the variational wave-function 
without both the long-range Jastrow factor and the pro- 
jector onto a fixed number of particles. 

In Fig. n we plot the variational parameter A„ar as a 
function of U for J = 0.2t, 5 = 0.16 and for different 
values of V. For J > V, Ayar starts finite at C/ = and 
increases with U. For V > J, Ayar = at small U, in 
agreement with the Hartree-Fock results. More remark- 
ably above a critical Uc a finite Ayar appears. Namely, 
the normal metal at > J turns into a superconductor 
by increasing the on-site repulsion. Both results can be 
explained within the Fermi liquid picture provided by the 
Gutzwiller approximation, according to which the effec- 
tive J* acting between the quasiparticles stays essentially 
unrenormalized when U increases, contrary to the effec- 
tive Vt, , which is substantially suppressed with respect to 
its bare value V. Therefore, as U increases for J > V, 
the quasiparticle bandwidth gets reduced, the attraction 
staying unrenormalized, so that the dimensionless cou- 
pling A* increases, hence Ayar- If J < ^, a normal metal 
is stable until > J.^ ~ J, after which superconductiv- 
ity gets in. In our numerical study we have found that the 
inclusion of the long range Jastrow factor Q consider- 
ably improves the simpler Gutzwiller wave function and 
allows larger values of Ayar- However, as shown in Fig. 
^3 the qualitative behavior of A^q^ vs. U is reproduced 
already by GA. 

Within the GA it is possible to study explicitly the 
competing influence of both J and V on superconductiv- 
ity. Let us consider the superconducting contributions to 
the uncorrelated expectation values for nearest neighbor 
sites i and j, 

(*scs|n.'^j|*i3cs)sc = 2A|c (5) 

{■^BCS\S^■SJ\■^BCS)SC = -^A|c, (6) 

where Asc = |(*bcs|c1^c]_^|^'bcs)| is the order pa- 
rameter of the uncorrelated wave function. In the case of 
Eq.llSJ, this term derives from configurations in which i 
and j are both singly occupied, both doubly occupied or 
one singly and the other doubly occupied, with weights 
(5^, (1 — S)'^ and 26{1 — S), respectively, where S is the 
doping. On the contrary, Eq.lQ has contribution only by 
configurations where both sites are singly occupied. In 
the limit of very large U, the configurations with doubly 
occupied sites are projected out, hence only the contri- 
butions from singly occupied sites survive in Eqs. © 
and (|SJ|. This implies that (jSJ gets a reduction factor S"^ 
relatively to ©, so that the actual condition for super- 
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conductivity at J7 ^ cxd reads approximately 



Since, in the same limit, the wave function renormal- 
ization Z ~ 2(5, we indeed recover the desired Fermi liq- 
uid result that interactions involving the charge density 
operators get renormalized down by a factor Z"^ with re- 
spect to those involving spin operators. The above dis- 
cussion also shows that not all the pairing mechanisms 
are equally equipped against on-site repulsion. Indeed a 
weak coupling d-wave superconductivity might be stabi- 
lized even by a negative y at J = 0: an explicit attraction 
between charges. However, for increasing J7, the effective 
strength of this attraction would decrease as Z"^ so that 
A, ^ ZXq, hence would go down. 

The behavior of the variational gap Ayar as shown in 
Fig. ^ suggests a crossover from weak to strong coupling 
superconductivity as U increases. This is manifest by 
comparing Fig. Ell with Fig. Eb. In the latter the varia- 
tional energy gap is displayed for several C/'s, while in the 
former we plot the true long range order parameter A ljiq 
in the correlated wave function. A^ao is estimated on 
a finite cluster through the pair-pair correlation function 
/ as follows: 



A 



1 



LRO 



where 



y<ci 

X, 



tS±1,- 



■IT y±l,-a' y, 



>, 



(7) 



being evaluated near to the maximum distance \x — y\ 
available on a given size. Notice that / includes normal 
contributions fnorm, which nevertheless vanish in the in- 
finite size limit. In order to improve the quality of any 
finite size analysis, one should estimate fnorm to get a 
meaningful value of the true long range order parame- 
ter. We decided to approximate fnorm by the value of / 



calculated with the optimized wave function having the 
same form with the variational parameter Ayar equal 
to zero (see inset FigEK)- After this subtraction, size ef- 
fects are acceptable, at least for our qualitative analysis 
(see FigUi). 

The crossover region where the gap Ayar rapidly moves 
from small BCS-like values to much larger values corre- 
sponds to a maximum of the true order parameter, as one 
would expect in the intermediate region between weak to 
strong coupling superconductivity. The notably differ- 
ence with the latter is that in our model the crossover 
does not occur by varying the bare attraction A, but by 
increasing the repulsion U. 

The different behavior of the variational gap with re- 
spect to the true order parameter, which has been asso- 
ciated with the behavior of the pseudogap versus Tc in 
the cuprateslslO], has a clear explanation within the GA, 
where A^no is suppressed by the factor Z with respect 
to the uncorrelated Asc- Indeed, as shown in Fig. EJi, 
the quasiparticle residue Z, defined as the jump in the 
momentum distribution function along the nodal direc- 
tions, is a decreasing function of U tending to Z ^ 25 as 
U -*oo. 

However, as shown in Fig. ^jp, Z is not the reduction 
factor of the full quasiparticle bandwidth, which gets con- 
tributions also from J and V. Again, this is an obvi- 
ous result in the GA where the Hartree-Fock decoupling 
of the nearest neighbor interactions effectively generate 
hopping terms. In spite of that, the charge current vertex 
is still determined by the true hopping i, hence gets sup- 
pressed by a factor Z ~ 2S. On the contrary, spin current 
vertex does include a contribution from J and survives 
against the strong wave function renormalization Z. 

In conclusion, we have shown that strong short range 
correlations enhance or suppress pairing correlations if 
they primarily involve spin or charge degrees of freedom, 
respectively. This behavior is manifest at strong U, in 
agreement with slave boson approaches'!^ and numer- 
ical calculations!^ 13. but starts to appear already 
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FIG. 2: Superconducting order parameter I^lro a), and the variational gap A„ar b), as a function of (7 a.t V = Q, J /t = 0.2, 
5 = 0.16. The inset in a) shows the long distance pairing correlations for the non-superconducting state (see text). 




FIG. 3: Panel a): wave-function renormalization Z as calculated through the jump in the momentum distribution along the 
nodal direction. Finite size scaling from 50 to 162 is used to evaluate the jump Z in the thermodynamic limit for the VMC. 
Panel b): quasiparticle bandwidth normalized to its uncorrelated value. The VMC refers to the 50 site cluster, as finite size 
effects are small. V = G and J /t — 0.2 for both figures. 



at weak coupling. Indeed, a recent calculation within 
the Random Phase Approximation shows that the d- 
wave superconducting phase of model 1^ at ^ = gains 
more exchange-correlation energy than a normal metal, 
thus supporting the results here obtained by variational 
Monte Carlo. 

JVote added. When this work was completed we be- 
came aware of a preprint by F.C. Zhang 13] which con- 
siders the Hamiltonian ^ with V = within the GA, in 
the context of the gossamer superconductivity scenario 
recently proposed by R. Laughlin Q . The results of 
Ref. [iJl qualitatively agrees with our VMC data. 
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